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Magnetization manipulation is essential for basic research and applications. A fundamental question is, how 
fast can the magnetization be reversed in nanoscale magnetic storage media. When subject to an ultrafast 
laser pulse, the speed of the magnetization dynamics depends on the nature of the energy transfer pathway. 
The order of the spin system can be effectively influenced through spin-flip processes mediated by hot 
electrons. It has been predicted that as electrons drive spins into the regime close to almost total 
demagnetization, characterized by a loss of ferromagnetic correlations near criticality, a second slower 
demagnetization process takes place after the initial fast drop of magnetization. By studying FePt, we 
unravel the fundamental role of the electronic structure. As the ferromagnet Fe becomes more noble in the 
FePt compound, the electronic structure is changed and the density of states around the Fermi level is 
reduced, thereby driving the spin correlations into the limit of critical fluctuations. We demonstrate the 
impact of the electrons and the ferromagnetic interactions, which allows a general insight into the 
mechanisms of spin dynamics when the ferromagnetic state is highly excited, and identifies possible 
recording speed limits in heat-assisted magnetization reversal. 

The FePt Ll 0 alloy represents the most important material for novel concepts in magnetic recording due to its 
high magnetic anisotropy, which ensures long-time thermal stability of nanometer sized bits 1 . Thin films of 
FePt with perpendicular anisotropy and small grain sizes are the most promising candidates for heat-assisted 
magnetic recording, which could reach storage densities beyond 1 Tb/inch 2 . Patterning continuous FePt into 
individual bits 2 can in principle extend recording densities to 100 Tb/inch 2 . The ultimate magnetic recording 
applications will also require faster bit switching and a deeper insight into the processes involved. However, non- 
deterministic fractioning in ultrafast magnetization reversal can limit the switching speed in recording schemes, 
and thus has inspired fundamental research for more than a decade 3 ' 4 . Recently, a new concept of ultrafast all- 
optical magnetic recording with an unprecedented switching timescale below 1 ps was suggested 5 . This opened 
new possibilities to reduce the speed limit established by the spin-orbit coupling timescale to that governed by the 
much stronger exchange interaction. Here we show that in FePt the intrinsic dynamics surprisingly is connected 
with the electronic structure: the alloying of Pt, because of its low lying d-bands, reduces the density of states 
significantly 6,7 and increases the electron temperature significantly above the Curie temperature. We identify the 
large electron energies of the heated electron system in FePt Ll 0 , arising from its reduced density of states at 
around the Fermi level, as its origin driving the spin fluctuations 8 ' 9 . 

Even though CoPt 3 was among the first thin film systems investigated 10 since the discovery of ultrafast 
demagnetization in 1996 by Beaurepaire et al. 11 , most investigations were centered on samples with in-plane 
anisotropy. So little is known about the behavior of materials with high perpendicular anisotropy 12 " 14 . A notable 
exception of materials with perpendicular anisotropy is the ferrimagnetic CoFeGd, which was studied in all- 
optical ultrafast switching triggered by a single laser pulse 515 . The modeling of this mechanism involves the 
exchange interaction of the two spin subsystems 16 . Different models have been proposed 17 " 21 . On the other hand, 
recent work of Koopmans et al. 22 suggests the classification of materials as "fast" (or type I) and "slow" (or type II) 
based on the ratio T c l\i at , where \i at is the magnetic momentum per atom and T c is the Curie temperature. In both 
cases, there is an initial sub-picosecond fast demagnetization. However, in the first case the fast femtosecond 
demagnetization is followed by a magnetization recovery (as in Ni 23 ), while in the latter a second slower demag- 
netization takes place. The recovery occurs on the timescale on the order of 50 ps and more (as in Gd 24 ). 
According to this classification FePt should be regarded as a fast magnetic material. However, more recently it 
has been shown that in Ni both behaviors can be observed 25 . Thus, the question of whether the high anisotropy 
material FePt can behave as "fast" or "slow" under specific laser excitation is an open question. In addition, for 
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thin films and granular media, the contributions of spin currents to 
the ultrafast demagnetization dynamics cannot be neglected 26-28 that 
can be directly detected by the emission of THz radiation 29 . We use 
insulating substrates and cap layers to minimize these effects. A 
transition from type I to type II demagnetization is found, triggered 
by the laser fluence. Theoretically this is analyzed by mapping the 
properties of the electron system onto a Heisenberg-like spin system 
in a thermal ensemble that determines the timescales of the ultrafast 
fluctuations 30 . A coupling parameter /L, determined by the micro- 
scopic spin-flip transitions, couples them 31 . Via two scenarios, we 
calculate a different response of the electronic system to the laser 
pulse. We can pinpoint that the quasi- equilibrium values of the 
electron temperature above the Curie temperature is determining 
the response of the spin system. This is leading to the critical fluctua- 
tions and a slowing down, a fundamental limit to recording speeds in 
heat-assisted reversal on the ultrafast timescales. To enable progress 
in high-speed and high-capacity magnetic storage devices, a fun- 
damental understanding of these dynamic processes is required in 
the future. 

Results 

We have studied a 3 nm-thick continuous FePt thin layer and a 
7 nm-thick AgCuFePt-C granular recording media, shown in 
Fig. 1 a) and b). In the granular media carbon intercalates the mag- 
netic 5-8 nm grains to separate them magnetically 32,33 . In Fig. 1 c) the 
time-resolved magneto -optical Kerr effect (A0 K /0 KyO ) is shown for 
the granular FePt recording medium as well as the data for the 
continuous FePt layer for a small demagnetization of below 3-7 % 
M s , which was set to +/— M s beforehand. Besides the absolute scale, 
both can be described using identical parameters of the analytical 
solution of a rate equation model, shown as a continuous line: the 
microscopic mechanisms on the nanometer length scale dominate 
the dynamics on a femtosecond timescale. This leads to identical 
spectra for the continuous film and granular media. Moreover for 
the case of the recording medium having granular structure, the 
larger H s and different K u does not alter the magnetization dynamics. 
The energy scale of the magnetic anisotropy K u < 1 meV is too small 
to affect the dynamics on the timescale related to the energy scale of 
the exchange interaction. Simultaneously to the magneto -optical 
Kerr rotation of the probe beam, the time-resolved reflectivity is 
determined. The decay of the reflectivity signal is fitted to an expo- 
nential function before characteristic stress waves set in (Fig. Id)), 
according to R(t) ~ exp ( — t/T E ). The results are presented in Fig. 1 f) 
showing the evolution of the characteristic timescale t e that 
represents the relaxation time of the electron temperature. Experi- 
mentally, the value for t e ~ 1 ps is a typical value expected for 
transition metals and consistent with previous results 35 . While for 
oxides discrepancies of the Kerr rotation and the genuine magnet- 
ization can be found even on long timescales 36 , this is meanwhile 
accepted not to be the case for metals for time scales longer than the 
electron equilibration time (about 50-100 fs). 

For the continuous FePt layer, using moderate B -fields of 200 mT 
saturates the magnetization and we can extract the microscopic para- 
meters of the ultrafast magnetization dynamics of FePt in polar 
geometry for a strong demagnetization. In the experiment, the flu- 
ence of the pump beam is varied from 5 to 35 mj/cm 2 . The Kerr 
rotation is extracted from that of opposite external field direction 34 ' 37 
to get the absolute degree of demagnetization, presented in Fig. 2. 
The first rapid demagnetization, t m , occurs at a timescale of 0.15- 
0.38 ps, increasing with the laser pump fluence. A "kink" is observed 
here on the short timescales at zero delay known from preceding 
experiments, for example by Bigot et al. 38 . The remagnetization slows 
down as a function of the incident fluence, until finally a second, 
slower demagnetization with the absence of recovery is found above a 
fluence of 30 mj/cm 2 . Thus, a transition from type I to type II as 
classified in Ref. 22 is observed for FePt. 
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Figure 1 | Sample characteristics: schematics and sample structure of the 
granular FePt recording media (a) and thin film (b) as measured by 
transmission electron microscopy (TEM). (c) Ultrafast magnetization 
dynamics for both cases after femtosecond laser excitation for a small 
demagnetization. Solid line: analytical three temperature model to obtain 
t m - Both curves can be described with sets of identical parameters, (d) The 
reflectivity dynamics from which the exponential decay z E and (e) the 
electron temperature T e are obtained (laser fluence 5 mj/cm 2 ). (f) 
Relaxation time t e for the electron temperature and t m for the ultrafast 
demagnetization is given as function of the pump fluence. The dashed line 
marks the linear increase. 



To understand the behavior, we model the ultrafast magnetization 
dynamics under external laser excitation by a thermal process of 
electronic origin 23 . The model considers that within a timescale of 
the order of 10 fs the electrons are thermalized and can be described 
by a quasi-equilibrium electron temperature T = T e (t), coupled to 
the spin system via the the fluctuation- dissipation theorem. We 
implemented the modeling within a stochastic micromagnetic 
Landau-Lifshitz-Bloch (LLB) model. The LLB equation for the spin 
S represents a thermodynamic extension of the Landau- Lifshitz- 
Gilbert (LLG) equation and has been derived by modeling distri- 
bution functions for thermal ensembles 39 . The thermal bath is the 
electron temperature T = T e (t), and the coupling- to-the bath para- 
meter X is related to microscopic scattering mechanisms. For the spin 
quantum number S=l/2 this versatile equation is equivalent to the 
microscopic three temperature model (M3TM) proposed for the 
description of ultrafast processes 35 . Using the LLB with spin S = 
3/2 it has been shown by comparison to an atomistic model for 
FePt that this describes its thermal equilibrium properties, e.g. damp- 
ing and magnetization versus temperature, as well as the ultrafast 
magnetization dynamics 30 ' 40 . Note that in the classical case, e.g. for 
S — > oo, the LLB equation can be derived as a coarse-grained approach 
from a classical Fokker- Planck equation for a thermal ensemble of 
spins, whose individual motion follows the Landau-Lifshitz-Gilbert 
equation, coupled to a thermal source via stochastic fluctuations in a 
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Figure 2 | Ultrafast demagnetization dynamics of FePt: spin dynamics 

measured by the Kerr set-up at increasing laser fluence in steps of 

5 mJ/cm 2 from 5 mJ/cm 2 (upper curve) to 35 mJ/cm 2 (lower curve). The 

detail on the femtosecond timescale is shown with expanded scale on the 
left side. For the highest fluence the demagnetization still slowly progesses 
after the first, fast demagnetization is observed (type II material). Data is 
normalized to the magnetization value before the pump pulse arrives 
(negative delay). 

similar way. The LLB's advantage, as compared to atomistic simula- 
tions, is that it is far less computationally intensive and thus 
represents a thermodynamically consistent extension of classical 
micromagnetism. Therefore, to extend its capabilities describing lat- 
eral inhomogenities (e.g. lateral thermal gradients or domains), here 
the LLB is extended to a thermal micromagnetic LLB. Combining 
individual thermal macrospins in a discrete array allows to model 
systems up to the micron- size. Thereby, each cell represents a ther- 
modynamic average over atomistic magnetic moments. Thus, ther- 
mal ultrafast heating is taken into account within each single cell by 
the longitudinal relaxation, equal to t m . However, this micromag- 
netic LLB allows describing locally a different thermal demagnetiza- 
tion, as well as a rotation against each other. In the model 900 thermal 
macrospins with a lateral cubic discretization of A = 3 nm with 
periodic boundary conditions are combined (Fig. 3). In detail, every 
single macrospin within the LLB formalism m, = M//M e (0) is 
described by the following formalism derived by Garanin 39 ' 41 ' 42 : 

M z is the magnetization (thermal average of atomistic spins over the 
volume V at temperature T), M e (T) is its thermal equilibrium 
value and M e (0) is the maximum at T = OK. The temperature 
dependence of M e (T) is evaluated via the Brillouin function with spin 
quantum number S = 3/2 43 . The effective field = H + H,-^ + 
H; 5 £x + H z j is comprised of applied, anisotropy, micromagnetic 
exchange, and internal exchange fields. The micromagnetic aniso- 
tropy field Hf a = — — (jni,x£x + Wi,y£y) is determined by the per- 
pendicular susceptibility x L = dm/dH±. Experimental 44 and 



theoretical 45 results report that in FePt the anisotropy K scales 
with magnetization as ~ m(T) 21 . Thus, we use x± = 
M 5 (0)/2iC(0)m(T) 01 . The micromagnetic exchange field is defined 

as 30 
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where; goes over neighboring elements and A,-(T) is the micromag- 
netic exchange stiffness. For FePt, we derived its temperature 
dependence from atomistic calculations using a FePt Hamiltonian 
with parameters taken from local spin density approximation 
(LSD A) 46 . It scales with the magnetization as ~ m(T) 176 . Here we 
use A i (T)=A(0)m(T)] J6 , where A(0) = 2.2 • 10" 6 erg/cm. The 
internal exchange field H zJ results from the thermal average of ato- 
mistic spins. At low temperatures, it is responsible for keeping the 
magnetization magnitude constant. At the same time this term is 
responsible for the critical behavior approaching ~ T c , which can 
be seen in the following expression: 
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In this expression a s = 2(S + 1) 2 /([S +1] 2 + S 2 ) and X\\ =dm/dH l 
represents the longitudinal susceptibility. The latter follows from 
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where B'Q stands for the derivative of the Brillouin function. The 
relationship between the internal exchange parameter / 0 (also related 
to A(T = 0 K), see Ref. 46) and T c is given by T c = S(S + l)J 0 /3k B 
where k B is the Boltzmann's constant. The stochastic fields & y± and 
Cud are given by 47 
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where i and j denote the macrospin number and k and / denote its 
Cartesian components x,y, and z. Finally, the longitudinal and trans- 
verse relaxation parameters are 42 



2T 2q s 
3T C sinh(2g 5 ) 



Si =X 



tanh(<^ 5 ) T 
q s 3Tc 



(7) 



where X is the microscopic relaxation parameter that couples the spin 
dynamics to the electron temperature T = T e (t), defined by the 
microscopic spin scattering rate /L, and q s = 3T c m e /[2(S + 1)T]. 

Via microscopic spin scattering rate k and the relaxation 
dynamics, the magnetization dynamics is coupled to the electron 
temperature T = T e (t) as described above. However in turn, the 
electron temperature within the two temperature model (2T) is 
coupled to the lattice temperature T ph (t) via rate equations, which 
determine its evolution with time t after the pump pulse arrived: 



C e — r^- = —G e -ph(T e — T p h) +P(t) — C e — 



dt 

Cph ^ = G e —ph (T e Tph) 



(8) 



Here C e and C ph denote the specific heat of the electrons and the 
lattice, respectively, G e _ ph is the coupling constant determining the 
energy exchange between the electron and lattice systems, and T ph = 
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Figure 3 | Density of states of Fe in (a) and FePt Ll 0 compound in (b): Pt alloying reduces the density of states at the Fermi level by lower lying 
d-bands of Pt 6 ' 7 . By blue, red and grey are Fe majority, Fe minority and Pt states are given respectively. In darker color the occupied states are given, 
(c) Illustration of the microscopic interaction of the heated electrons T e and the spin system via spin flips coupled by X for both cases, (d) Ultrafast 
dynamics for a hypothetical material having a large density of states and FePt using a small density of states (related to y e ) for a high laser fluence. In both 
cases the same reflectivity dynamics (z E ) and maximum demagnetization is calculated, but the resulting ultrafast magnetization dynamics is different, 
(e) Implementation of the micromagnetic Landau-Lifshitz-Bloch (LLB) model: the magnetization is described by the average over 900 thermal 
macrospins with a lateral cubic discretization of A = 3 nm with periodic boundary conditions. Each cell represents the thermodynamic average over 
atomistic magnetic moments (shown schematically on the right side). The thermal ultrafast heating is taken into account within each single cell by the 
longitudinal relaxation M(t) and the individual macrospin allows a rotation. The lateral discretization allows the implementation of more general lateral 
inhomogeneous excitation (e.g. temperature gradients). 



340 ps is the heat diffusion time to the substrate. For Q, the free 
electron approximation is used resulting in C e = y e T e , where y e is 
determined by the electron density of states. C P h is set constant, since 
FePt has a Debye temperature well below the room temperature 
T r0 om> The laser absorbed power is defined by P(t) = I 0 F 
exp [ — (t/Tp) 2 ] proportional to the laser pump fluence F. The time- 
resolved reflectivity reveals that the change of electron temperature 
depends linearly on the change in reflectivity 48 . Thus, we assume that 
the lacking parameters of the 2T model can be extracted from the 
experimental data. The thermal diffusion timescale to the substrate 
T ph is obtained from the long-term magnetization behavior. The 



proportionality constant I 0 is estimated by fitting of the experimental 
demagnetization value at 30 ps and the coupling-to-the bath para- 
meter X via matching the maximum demagnetization value. 
The microscopic parameter X is related to the spin-flip proba- 
bility. Assuming that the main spin-flip scattering mechanism 
is the same that is responsible for the energy dissipation pro- 
cesses on a longer timescale, X is connected to the macroscopic 
Gilbert damping. We used 0.01 to 0.1 as a starting point for 
FePt 30,49 " 51 . This value scales with the spin-orbit interaction 51 . 
But generally the values of the Gilbert damping are found to be 
larger than the electronic coupling-to-the bath parameter X. 
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Figure 4 | Simulation of the electron temperature T e , shown as a function 
of the laser pump fluence (from 5 mj/cm 2 (upper curve) to 40 mj/cm 2 
(lower curve), in steps of 5 mj/cm 2 . The 2T model is based on the set of 
parameters presented in Table I. The parameters are extracted from the 
reflectivity dynamics (Fig. 1). Within the shaded area marked in the left 
panel, the electron temperature exceeds the Curie temperature (dashed 
line). 



Additional contributions to the Gilbert damping are magnon- 
magnon scattering and inhomogeneous line broadening of the 
ferromagnetic resonance 34 ' 52 . 

The determination of the material specific constants such as y e and 
G e -ph is crucial for a proper simulation of the demagnetization pro- 
cess. Two approaches were performed, resulting in a high and a low 
electron temperature profile. To model the profile for the FePt, first 
G e _ p h = 1.5 • 10 17 W/m 3 K is assumed in the range of Cu, Mo, and 
Pt 53 . Then, the analysis of the reflectivity relaxation rate gives y e = 
110 J/m 3 K 2 . This value is of the order reported for Au and Cu 53 , but is 
much smaller than the corresponding value for Ni and hence pro- 
duces a high electron temperature 23 . In the second case of a hypo- 
thetical material, a coupling constant G e _ ph = 1.8 • 10 18 W/m 3 K, 
similar to Ni 23 , is assumed which gives y e = 1700 J/m 3 K 2 . This y e is 
the proper value for transition metals with a large density of states at 
the Fermi level due to crossing d-bands. As a consequence, a lower 
electron temperature is reached (Fig. 3d)). 



Discussion 

The results for the integration of the set of the LLB equations, 
coupled to the 2T model, with the set of parameters for FePt are 
presented in Fig. 5 for all fluences. As in the experiment, the simula- 
tions show a transition between type I and type II behavior. A fair 
agreement is found for the evolution of the demagnetization values 
AM/M(300 K) = 0.05-0.7. The demagnetization time t m are 
matched, and range in between t m = 0.2-0.35 ps and, within this 
fluence range, they show an increase (Fig. If)). Because of the large 
electron temperatures, a small k = 0.01 is sufficient to reach the 
demagnetization values. Contrary to this, in case of the integration 
of the LLB equation coupled to the 2T model with y e = 1700 J/m 3 K 2 
for the hypothetical material, the demagnetization time t m stays at a 
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Figure 5 | Ultrafast demagnetization dynamics obtained by integration 
of the LLB micromagnetic model coupled to the 2T model (Fig. 4). The 

model data are shown with increased laser fluence from top to bottom in 
steps of 5 mj/cm 2 from 5 mj/cm 2 (upper curve) to 40 mj/cm 2 (lower 
curve). The detailed view on the femtosecond timescale is shown with 
expanded scale on the left side. 

constant value and is always followed by the remagnetization within 
several picoseconds (see Supplementary Materials). Here a large k = 
0.1 is necessary to reach the experimental demagnetization values. 
The difference can be seen in Fig. 3 d) where both scenarios are 
presented directly opposed to each other. We obsereve that the sec- 
ond scenario does not produce a transition between type I and type II 
behavior. The origin of the high electron temperatures is the reduced 
density of states of FePt by Pt alloying. Because of lower lying d- 
bands in Pt as compared to Fe (Fig. 3 a, b), the ferromagnet Fe 
becomes more noble in the FePt Ll 0 compound 6 ' 7 . This is the under- 
lying reason for the electron specific heat being more similar to Au 
and Cu with a small y e . For the hypothetical scenario with large y e , 
only below 1 ps a T e above the Curie temperature T c is found. For the 
full calculation for all fluences refer to the Supplementary Material. 
We conclude that the transition from type I to type II demagnetiza- 
tion is defined by a critical temperature that the electrons have 
reached. The microscopic picture behind is that the electron system 
drives via spin-flips the magnetic system into a region, where the 
ferromagnetic correlations are strongly disturbed as schematically 
shown in Fig. 3 c). The role of the coupling parameter k is only to 
determine the rate of spin-flips 31 . Microscopic spin fluctuations are 
responsible for the decrease of the micromagnetic exchange stiffness 
A(T) and thermal behavior of M e (T) within the LLB model. The 
crucial effect, however, is that they slow down the magnetization 
response when the electron temperature is approaching the Curie 
temperature. In the LLB theory this is taken into account by the 
critical behavior of the longitudinal susceptibility, diverging at the 
Curie temperature. The effect is seen by looking at the shaded regions 
in Fig. 3 d) and Fig. 4. Only if the electron temperature T e stays near 
the Curie temperature T c for several picoseconds this transition in 
demagnetization behavior is found. This is a collective effect, consti- 
tuting the magnetitization dynamics: from the theory of phase tran- 
sitions we know that the dynamics is dominated by spin fluctuations 
and divergence of the correlation lengths. This leads to a slowing 
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down of correlation times, and connected to that, a longer demag- 
netization time t m 40 ' 54 . The transition to the type II behavior occurs 
when the timescale t m , here given in an extension valid approaching 
T c , exceeds the timescale t e 22 : 

Vat 



2ykk B {T c -T e 



>?E- 



(9) 



For the FePt case, a delicate balance between the driving electron 
temperature T e , determined by the small specific heat of the electron 
system for FePt, and the other parameters, drives the observation of 
the type II demagnetization in a small window of experimental para- 
meters. This is in contrast to Ni, where for the same parameters only 
a slowing down is observed 23 . 

However also certain interesting discrepancies are observed. In 
comparison to the model (Fig. 5), in the experiment (Fig. 2) the 
transition from type I to type II demagnetization takes place for a 
lower fluence. To investigate this further, we plot the demagnetiza- 
tion dynamics for the highest fluences, now normalized to the value 
of M s at room temperature when the pump -beam is blocked. This 
value can be determined by the hysteresis curves given in Fig. 6. The 
demagnetization relative to M s is lowered further with pump fluence 
and the two step behavior of the magnetization changes to a decrease 
within the first picosecond and stays on a plateau. This discrepancy 
between experiment and theory is related to the decrease of the 
magnetization at negative delay. The magnetization does not recover 
in between two laser pulses (e.g. 4 fis pulse repetition). These non- 
reversible effects are also mirrored in the strong reduction of the 
coercive fields, ideal for heat- assisted writing, meaning that the 
reversal process is significantly altered. We do not have a full under- 
standing at the moment. Two effects could be the origin that the 
magnetization dynamics slows down so strongly and stays at its 
value. The heat is not conducted away and the repeated laser pulses 
increase the basis temperature, or that the magnetization itself does 
not recover in between two pulses completely. But it will be most 
interesting to explore the origin and importance for a reliable heat 
assisted recording technology in future hard-disc drives. 
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Figure 6 | Ultrafast demagnetization dynamics experiment at highest 
fluences for 30 mJ/cm 2 to 40 mJ/cm 2 . Here the time-resolved data is 
normalized to its magnetization value with blocked pump beam. The 
hysteresis curves for 5 mJ/cm 2 and 40 mJ/cm 2 are shown for blocked 
pump-beam and negative delay of —5 ps. They demonstrate the 
irreversible parts of the laser induced dynamics present before the new laser 
pulse arrives (in a time scale of 4 /us pulse repetition period), becoming 
prominent for 40 mJ/cm 2 . For 0.5 ps and 2 ps delay the hysteresis is scaled 
to the maximum demagnetization, demonstrating the size of the 
demagnetization in total (reversible and irreversible contributions). 



Table I | Overview of parameter set I for the simulation in the case of 
high electron temperature for FePt and set II for a hypothetic mater- 
ial with a large density of states at the Fermi level, the hypothetic low 
electronic temperature case 

Parameter set y e X l 0 C p h T p /, S 

(J/m 3 K 2 ) (s" 1 ) (W/m 3 K) (J/m 3 K) (ps) (h) 



110 0.01 5.0- 10 1 ' 
1700 0.1 3.0-10 i; 



1.5 • 10 17 3.7- 10 5 340 3/2 
1.8- 10 18 3.3- 10 6 340 3/2 



In summary, by means of time- resolved Kerr magnetometry we 
have investigated ultrafast magnetization dynamics in FePt thin films 
with perpendicular anisotropy. Our results indicate that the electron 
temperatures reached by the laser heating play a crucial role in the 
character of the ultrafast demagnetization in FePt. The measure- 
ments reveal a transition from type I to type II behavior. Our experi- 
mental results are modeled in terms of the micromagnetic LLB 
model, coupled to the 2T model. We discuss two scenarios with 
different specific heat of the electron system. Within this framework, 
we find that a transition to type II behavior is a consequence of the 
high electron temperature. We identify that at large pump fluences 
the resulting electron temperature remains close to the Curie 
temperature and is leading to critical magnetization fluctuations 
responsible for this transition 40 . This non- deterministic spin 
dynamics is responsible for a speed limitation of the magnetic 
response to the laser pulse. Note that this is not only determined 
by the ultrafast optical pulse, but also by the nature of the FePt's 
density of states at the Fermi level, which defines the increase in 
electron temperature. From our experiments we conclude that a 
careful adjustment between electron heating is essential, while the 
anisotropy energy of the system does not play a crucial role for the 
response time. We identified the nature of the electronic density of 
states as a tuning parameter for the strength, a spin system will 
react to the laser power. Our results open possibilities for ultrafast 
control of the demagnetization in FePt, the most promising can- 
didate for future magnetic recording. Importantly, we have shown 
that we are able to manipulate the degree of demagnetization and 
its ultrafast rates. We propose that for efficient writing the degree 
of heating and its speed have to be balanced by varying the amount 
of the energy deposited. 

Methods 

Fabrication. The continuous FePt thin layer and AgCuFePt-C granular recording 
media have been prepared by sputtering. The structure from bottom to top is glass/ 
NiTa/MgO/FePt/Si0 2 and glass/NiTa/MgO/AgCuFePt-C/carbon overcoat 
respectively. The granular film has small amount of Ag to reduce the FePt Ll 0 ordering 
temperature, and a small amount of Cu to lower the T c . Carbon intercalates the 
magnetic 5-8 nm grains in a specific heating cycle 32 . The continuous film sample has 
a Curie temperature of 650 K, a saturation magnetization M s = 1070 emu/cm 3 , a 
maximum anisotropy constant K u>max = 2.2 • 10 7 erg/cm 3 , and an average anisotropy 
constant K u>av = 1.4 • 10 7 erg/cm 3 . Hystersis curves are given in the Supplementary 
Materials and reveal for the granular recording media [i 0 H s = 2.4 T and for the 
continuous FePt film |i 0 H s = 150 mT. 

Experimental set-up. In the all-optical pump probe experiment the magneto -optical 
Kerr rotation of the probe beam is measured 34 . The pulse has a full width at half 
maximum of x p = 40 fs and a wave length of X = 800 nm (Ti:Sapphire amplifier 
system RegA 9040, Coherent). The repetition rate of the laser system used was 
250 IcHz. The Kerr rotation is extracted from that of opposite external field 
direction 37 . To get the absolute degree of demagnetization, the Kerr signal is scaled to 
hysteresis measurements at two states of reference; one at negative delay (6 K ~ M z 0 ) 
and the other at a time delay that shows the lowest magnetization (A6 K , min ~ 
AM z>min ). 

Thermal modelling. To describe the ultrafast magnetization dynamics, we model a 
thermal process of electronic origin 23 . The model considers that within a timescale of 
the order of 10 fs the electrons are thermalized and can be described by a quasi- 
equilibrium electron temperature T e (t). A multi-macrospin model is used with cubic 
discretization elements with a lateral size of A = 3 nm (volume V = A 3 ). Within each 
cell the Landau- Lifshitz-Bloch (LLB) micromagnetic formalism is applied 41 . For the 
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present simulation, a system of 30 X 30 X 1 macrospins with periodic boundary 
conditions in x and y directions is used. 
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